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AEFREITATIONS  A^H  S^KEOLS  TO  12  USED  IN  TiUS  KETORT 

fR  Wiiiii  aspect  ratio 

Cd  Aircraft  total  drag  coefficient 

Cq.  Induced  drag  coefficient  •■ 

Crv  Drag  coefficient  at  sero  lift 

^o  • 

C]»  Thrust  specific  fuel  consumption 

CF        .    Centrifugal  force 

D  Total  aircraft  drag 

d  Total  differential  operator 

Partial  differential  operator 


c^ 


D^  Induced  drag 

Dq  Drag  at  zero  lift 

e  Airplane  efficiency  factor 

F  Total  effective  thrust 

F*  Thrust  due  to  momentum  change  of  air 

f  Effective  thrust  available 

g  Gravitational  acceleration  constant   (32.17  ft/sec^) 

h  Altitude 

L  Aircraft  lift  > 

m  Mass  of  air  and  fuel 

K  Mach  number 

q  Dynamic  pressure  (q  ■  1  p  V^) 

S  Wing  surface  area 

t  Time 

u  Engine  exhaust  veiocitj  relative  to  the  origin 

u  Exhaust  velocity  relative  to  the  Lody  (V  -  u) 


V  Velocity 

Vj,  Speed  of  sound 

V  Velocity  derivative  with  respect  to  time 
W  Total  weight  of  bod;,  or  aircraft 

2W2  Weight  of  fuel  expanded  from  point  1  to  2 

X  Horizontal  distance 

'i  Flight  path  angle  relative  to  the  x-axis 

^'  Flight  path  anr-le  derivative  with  respect  to  time 

yO  Air  densitj' 

ui)  Angular  velocity  in  pitch 

Cl  Aircraft  total  lift  coefficient 


In  this  rerort  the  r^rincin-il."?  of  varl£:t,5.onnl  crlruir.s  .v.'e 
ap.iii',-'.  i,;   .;.:rxv  .qaatijoG  for  a  fli^jiit  path  ropi'ci;>_;.t,.i\    either 

naxl-iun  range  for  a  :;iv:;i  fuel  expenditure  of  iTu.aun..'j-  fuel  e.c  a-e 

for  5    ^i-ren  ran,:e.     The  ecj-iations  ro  deterrdnsd  nnrr.it   rolnh'^n  of 

,e  proLien  witiiout  nakinj  tl-ie  norr-^.l  r.'^ctrictive  assun;)tionr;  o:  lift 
equals  weight  and  thrust  equals  drag,     Tiiis  ocrmits  analytical  treat- 
ment of  flight  attituJeL  ran-:in--  f-ror.  i^vel  fli-;'-t  to  vt-rtir.-ii  fli  ■-'>■{. 
and  scceieiT-ulonG  up  to  Etruct'Ui'^i  or  pilot  iir;;iti.v:-o:;yo     v;ita  zach 
a   general  trer.tnent  zf  aircraft  motion,  it  will  Lr.  possille  to  conpute 
.rch.^dulor  for  naxir.uii  v^n:-Q  fll  -'-t  o^tlis   corx.encin  ;  ct  en-    ^^'r  5  *  (--Vj - 
i; -i-ai-iaLl'i     co:TiLin:;.tion  of  aioj  :.;  -,   c  jo^^d,  and  cl:.3tanc3. 

The  equatioai  oittained  for  these  extrer^l  fli;_;ht  pathc  are 
fiTt  rn-i  rrcon.'l  order  nonj-ine-ir  diff^r-^ntial  equatioir-.     N'-'-'-'ricai 
soL'i:.-  -.-,-    oi    taj.s  system  are  shown  to  Lc  pocsille  usin^j  etep'.risc  inte- 
gration tecliniquese 

Clar-Gxcal  'jaiMnetorr   ^ot     A'.-'r^i'irr'  rriri^e  f  1  ^  '■:''  •■^''-"  ■r'.-.'r'ifi-" 
derived  froiv.  L.;e  pat;i  equaticnc  Lj   cpaci.  ^in,-,'  level  fil^  .;  , 

indicating  tiiat  ciassicr.i  i'an_-'--^  equations  are  restrictive  apolicr.tionr 
of  the  nore    •■'ncr-il  equp,t' n-'c;   ■^r^r'i-.^r-'.      ■  tills  thesis. 

T:a.:j   investir^ation  ws   conducted  in  partial  fuifuilnent  of 
the  requirements  for  the  degree  of  ^aster  of  Science  in  Engineering 
hi  Lieutenant  Corinandar  Ch?.ri'-r  v„  vv-itr:,  usil,  and  Lieutenant  Joh:!  i;,  Tlorii^:,, 
USM,   d-.ur'ing  the  VJCU  -  1955  ac^ideriic  :.ear  at  the  Forrectal  Rec^irch 
Center,  Princeton  Univvarsit: ,  Princeton,  IJg;;  Jersey. 


The  prollem  of  dct'.riT:.lnin^  tne  raximun  range  schedule  of 
an  aircraft  with  n   s-^'c'fied     '"  el  lor.-]  hr.s  rocrAved  ver?    little 
attenti.on  since  the  Ire^^iet  dcy.-iopment£  r.in:  tne  r.arc  recent  v?ork 
on  cruise  clin.Ls  for  j-.l  t;,  ;■=;  aircraft.     TuGse  r.cthodc  are  vcr;   rectric- 
■f  •  ,-r.  -i  n   •  K-.  ^   the;^  appl:    principally    .      _^v- rrronujnc^-  si^cr.ift.     .*n: 
inrn'overfient,  eitner  in  aircrcft  dosij?!  or  in  fligiit  patn  scredui^c,  v/liich 
would  -Inrrease  the  ran^^e  for  a   ,^iven  ar.ount  of  fuel  burnsd  in  f.n  rir- 
cr?*;-  ^^t  enthusiapMc'-ll:'  recf^'w'  b'    -"Gci-ners,  tacticisnc, 

and  r.n;  one  concerned  with  aviation,     Iji  tir:  follo-.'ia^  pa^^r-^s  a  nethod 
fo-*^  .-v^t.  nnination  of  flijht  pat:    schedule::;  for  maximun  ran^e  will 

oo"'l  vr'-iic'Tfj   t!ie    -  '•!; '^'^"^  ■    ^*  ;'>1,   can  nr'j/'uce  ci'niific?  nt  inorove- 
'  'i.'i.ts  over-  an^'  method  now  in  use-. 

Kan's  intuitive  sense  x^ill  tell  hira  that  for  alj   po^-.Ule 
pc-.t'7.c    Ln^aoe,  t'\ere  vjill  be   ■■^""■^  -[^orticUiar  path  vihich  vrill  r-^^-^y  'rrn 
furthest  for  a  ;-;i'/en  fuel  exj^eaditure.     This  ct-..':r\2  sens=  tells  us  thfit 
for  r      i  vcn  dict-rnce  traveled  onlj-  one  path  trill  yield  the  iminijnuin 
fuel  consurqption,     Restrictin-;  tlie   :;enc:r£il  case  + -'  the  o-^-r.-rin-  diri^n- 
sions  and  capabilities  of  modern  aircraft,  it  can  le  seen  thrl  for  every 
flight  path  of  vTtriable  altitude  and  speed,  there  will  le  one  parti- 
cvilar  schedule  of  altitude  ar  '  <?  -=^..-:  ...rPiicn  v/ill  :>ve  the   greatest 
range  for  a  given  fuel  consiu:ption.     There  vfili  also  !  e  one  flight 
schedule  ■•■/Inich  vrilll  give  the  rrdninum  fuel  coasuription  for  a   '^ven 
rr-'    ,-o     I  i^o  1  illustrates  t^.is   ■-—^'■Icr-.Liio     The  nost  (general  air- 
r-^aft  problem  will  ha/c  loo3?,lj    detined  initial  conditions.     The  path 
"or  rninJTiun  fuel  co.;ld  le  started  fvom  ap.\   cor.ibination  of  inital  speed 
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ci  isfon*:.^ 


Fij.  1 


F;j.  ^ 


and  altitude,  such  as  pt.  A  on  fi^.  1,  and  fiovm  a  given  distance  x^. 
The  schedule  of  altitude  and  speed  would  be  detomined  by  the  inter- 
rolations  of  the  \'ariabloE  of  the  proLleLi  bounded  by  the  capabiiitisE   of 
the  particular  aircraft.     The  pVoblen  could  just  as  well  be  started 
fonn  the  instant  of  take-off,   say  pt,  E,  and  thus  cover  the  entire  flight 
from  the  take-off  until  over  3r]_, 

Cora.iencing  a  maxLmun  rf  ngo  or  rciniriun  fuel  sche<^ule  from 
take-off  Is  a  conpletely  new  concept,     Clas3ic?.l  range  solutions  do 
not  attempt  to  solve  for  an  optirium  flight  path  from  point  of  take-off 
to  optimun  cruising  altitude,  but  rather  specify,  for  turbo-jets, 
£  r-nebulous  "njniinuin  fuel  climb"   to  the  highest  altitude  r^t  vrhich 
speed  for  (<^LVCD)niax  ^^  equal  to  V^^^^.     The  aircraft  then  cruises 
at  that  altitude,  slowly  climbing  as  the  weight  is  decreased,  thus 
increasin-  the  altitude  at  v/hich  the  speed  for  maximiim  Cy^  /Cr)  is 
equal  to  the  naximuiii  speed. 

This  "minimun  fuel  climb"  will,  if  properly  determined, 
achieve  a  certain  altitude  with  a  minimum  fuel  expenditure,  but  Uiis 

is  not  the  primary  concern.     The  problem  is  to  travel  a  distance  xi 

f 

with  minimum  fuel  consumed.     If  the  variables  of  the  problem  call  for 

a  climb  to  a  specified  altitude,  then  that  cliirib  should  be  executed  to 
give  maximum  distance  for  a   given  fuel  expenditure  to  get  to  that 
altitude.     As  long  as  the  flight  nath  stays  within  the  limitations  of 
the  aircraft,  however,  there  should  be  a  single  path  with  no  dis- 
continuities for  minimum  fuel  path  for  a  given  distance.     Thus 
looking  at  fig,   2,  it  can  be  seen  t'lat  to  travel  from  pt.  a  to  some 
pt,  Si  a  distance  Xj^  avjay  there  must  be  some  minimum  path  ag  which  uses 


i 
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less  fuel  than  the  path  i^]?n  wliich  is  corrcposed  of  two  parts,  even 
though  ^  aad  Ic  aro  separate  patlis. 

Classical  range  equations,  moreover,  are  not  generally 
applical le  or  accurate  for  present  da^  high  performance  aircraft. 
Classical  ran^je  performance  is  bc.sed  on  assuming  that  the  flight 
path  angle  (l()  is  snail  such  that  its  cosine  is  considered  one, 
(this  assunes  that  L  ■  W)  and  that  flight  path  accelerations  along 
the  path  vjere  ignored  (i.e.  Thrust  ■  ''"^rag).     From  a  review  of  the 
literature  on  the  culject,  it  was  determined  that  ever;/thing  written 
to  date  on  the  range  of  aircraft  has  been  based  on  these  simplifying 
but  verv  restrictive  assuirpti  ons.     With  modern  high  performance  air- 
craft, flight  path  angles  up  to  and  including  90  degrees  are  possible 
and  flight  path  accelerations  can  easilj^  exceed  1,0  g.     In  fact,   in 
considering  the  gener.'il  problonof  powered  flip;ht  (rockets,  missels, 
VTO,  etc.),  the  convenient  and  conventional,  condition  of  steady  level 
flight  is  almost  the  exception  instead  of  the  rule. 

It  is  the  intent  of  this  work  to  indicate  methods  lased  on 
the  principals  of  variational  calculus  for  the  detsnr.ination  of  schedules 
for  minimum  fuel  or  maximum  range  flights  for  all  tjpes  of  powered 
flights.     It  is  intended,  more  specif icallv,  to  develop  equations 
for  the  mininum  fuel  paths  of  turbo-iet  t:,pe  aircraft.     This  differs 
from  bh-'  ciasGic  problem  in  that  the  schedule  begins  at  take-off  and 
minir.izes  the  fuel  expenditure  li'oir.  tiiat  point  to  a  distance  Xq_,  taking 
int'j  consideration  possible  large  values  of  flight  path  angle  and 
ronsider-ible  flight  path  accel^-rrtions.     In  other    words  ,  the  oroblen 
considered  nere  is  the  broad  general  problem.    It  uses  accurate, 
realistic  conditions  and  minimizes  fuel  consuirption  for  the  entire 


flight  path  from  "wheels  up"  to  "over  tar^^et". 

It  is  hrjrdly  necessary  to  discuss  the  iriportance  of  reduced 
fuel  consuiTDptian  since  its  advantages  are  so  readilj  apparent,  but 
to  help  keep  these  advantages  in  mind,  a  few  of  then  will  be  listed: 

(1)  Increased  useful  payload  (bombs,  cargo,  passengers,  etc) 

(2)  Greater  range  or  combat  radius  with  ttie  same-  paj'load 

(3)  Koro  economical  operation  overall 

(k)   Complete  specification  of  a  flight  schedule  which  could 
be  incorporated  into  an  automatic  control  for  a  manned 
or  pilotless  aircraft. 
The  adaptability  of  these  and  other  advantages  to  all  phases  of 
aviation  is  also  readily  apparent.  The  long  range  bomber,  the  short 
range  fighter,  the  commercial  cargo  and  passenger  carrier,  the  pilotless 
aircraft,  and  practically  all  other  vehicles  traveling  tiirough  the 
air  could  adopt  flight  schedules  determined  in  this  manner  to  increase 
their  effectiveness  and  their  efficiency, 

A  specific  application  for  which  this  work  could  be  im^. 
mediately  useful  is  the  determination  of  entry  and  escape  paths  for  a 
supersonic  bombing  mission,.  Those  parts  of  the  entry  and  escape 
path  not  determined  hi   other  tactical  considerations  would  most  prob- 
ably be  flown  at  full  power  including  afterburner  and  following  a 
path  to  give  best  fuel  economy  at  maximum  power.  The  bounded  thrust 
path  equations  deveiqjed  later  in  this  thesis  apply  directly  to  this 
problem. 

In  this  and  many  other  applications,  Kach  effects  must  be 
considered.  For  sinplicity,  however,  Mach  effects  have  been  omitted 
from  the  deri/ations  that  follow,  but  have  been  treated  separately  in 
Appendix  B» 


PART  TWO:   D2VJL0B  E^FT  CF  S^^^ATIONS 

The  range  of  an  aircraft  is,  as  we  have  hjpothesised,  a 
function  of  its  filrht  path.  Its  flight  path,  in  turn,  dqjends  on 
vp.rialles  such  as  amount  of  fuel,  velocity  aiaig  path,  flif^ht  path 
an-jie,  accelerations  along  and  normal  to  the  path,  and  thrust  spec- 
ific fuel  consumption.  These  variables  are,  in  turn,  functions  of 
each  other,  drag,  tlirust,  lift,  weight,  and  air  density.  Relating 
these  variables  among  sufficient  eq^oations  to  solve  for  a  minimum 
fuel  path  is  the  problem  of  thif  thesis. 

A)  EQUATIONS  OF  NOTION 

The  equations  of  motion  pertinent  to  the  problem  are  the 
force  balance  equations  of  a  body  moving  in  a  vertical  plane.  This 
motion  will  include  accelerations- along  and  noiTmal  to  the  flight  path, 
as  in  maneuvering  flight,  since  these  accelerations  will,  a  priori, 
affect  the  loads  and  thrust  required,  and  therefore  affect  the  range 
(or  fuel  consurption).  The  forces  acting  on  a  body  moving  along  an 
arbitrary  path  in  tvo  dimensional  space  are  shown  in  fig.  3» 

If  this  body  moving  i^  space  is  a  conventional  aircraft, 
the  forces  Fi,  F2,  and  Fo  can  be  expressed  asi 

Fi  The  total  effective  thrust.  (F) 

F2  The  lift  of  the  moving  body.  (L) 

F3  The  total  drag  force.  (Body  drag,  D,  and  inertia,  m  dV  ) 
By  equating  forces  parallel  and  perpendicular  to  the  flight  path  to 
zero,  the  following  equations  are  obtained: 

(1)  Forces  parallel  to  flight  path 

2:  F„  -  0  -  Fi  .  F3  -  WsinV  (1) 


A 


h     1 

-1 


^ 


PisTAHce;     r  X 


F.g      3 


FoRces         ACTit^^      o  t*      Boi>Y 


; 


or  for  an  aircraf tj 

F*-D  -  d  (mV)  -  Wsin^  ♦  u  ^  .  0  (2) 

where  F,  D,  W,  T^  ,  u,  and  in  were  defined  previously  andt 

d(inV)  -  in  dV  ♦  V  An 
dt  dt  dt 

But 

dm  ■  -OpF  /j|N 

dt  g 

since  fuel  is  the  only  variable  mass  involved* 

Also     dV  2  V 
dt 

Now 

Y^Y       -  F*-D  -  mV  -»•  VGpF  -  Wsini^  -  uCxF  -  0  (U) 


Considering  the  definition  of  thrust,  it  can  be  seen  that  the  total 

effective  thrust  acting  on  the  aircraft  would  be  equal  to  that  inv- 
parted  by  the  air  through  the  engine  as  it  experiences  a  change  of 
momcntiim  plus  that  inparted  by  the  fuel  mass  leaving  in  the  exhaust 
due  to  its  momentum  change.  The,  thrust  due  to  the  leaving  fuel  is  a 
function  of  the  difference  between  aircraft  velocity  and  exhaust 
velocity  and  is  e^^ressed  by: 

(u  -  V)  Cj  F 

g 


Thus 


F  »  F*  -  (u  -  V)  Cj  F 

»  F*  •»•  Ue  Ct  F  (5) 

g 


The  equation  now  becomes 

F  -  D  -  W  V-  Wsin^  -  0  (6) 

g 

(2)     Forces  perpendiciilar  to  the  flight  path 

Jf^     -  0  »  Fg  -  CF  -  W  cosi  (7) 

or  for  an  aircraft 

L  -  CF  -  VJcosi    -  0  (8) 


But 


and 


then 


C  F  -  m  V  J 


u>    - 

dt 

HV- 

C  F 

-  W 

i 

v^ 

And  therefore 


(9) 


EF^  -  L  -  W  V^  -  Wcos/  -  0  (10) 

g 

Equations  (6)  and  (10)  above  are  the  general  eqijations  of 
motion  pertinent  to  flight  in  the  vertical  plane.  These  will  be  modi- 
fied later  for  use  in  the  range  problem  to  the  extent  that  derivatives 
with  respect  to  time  will  be  replaced  by  parametric  derivatives  with 
respect  to  either  weight  or  distance. 

B)  EQUATIONS  OF  MOTION  FOR  RANGE 

Before  attempting  to  maximize  the  general  equations  of  motion 
we  will  further  modify  them  to  be  more  applicable  to  the  problem  of, 

(1)  Maximizing  the  range,  x,  for  a  given  fuel  ejqjenditure,  2W2»  °^ 

(2)  Minimizing  the  fuel  e^q^enditure,  ]^W2for  a  given  range,  x* 


/\. 
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7or  maxi^^izing  rancyi,  x,   it  is  convenient  to  clianfje  variables 
as   follows: 

Define  d(    )/dJ         (    )» 

Then  V  -  J/  -  JV  ^  -    y»    dv 

dt       d.:  dt  ^ 

But  d7  -  -Crp  F 

dt 

Thergfore  V     -     -7«    C^  F  (^1) 

Also 

dt  d"     dt  dt 

?  -     -i<«      C<j.  F  (12) 

Substitutin-  (11)  and  (12)  into  (6)  an ^   (10)  r^ves: 

L  +  W  V  lf»   Ct  ?  -  :Jcos)(  -  0  (13a) 

i 

F  -  D  ♦  W  V    Or  F  -  :rsin^  »  0  (13b) 

The  integral  equation  ir^latinp;    x     and  VJ     in  simplest  form  is: 

Ac  d' 

Now  d:c  -  dx  dt 

.■?;•        dt  r^: 

But  *£  •  Vcos'i 

dt 
And  fron  befoi^ 

dt  -  -     1 

d:  Oy  F 


so 


^  »  -  Vcos 8 


Ct  F 


-  Vcos     dv  (1E;> 

Ct  F 


Frora  (l3b) 

F  «     D  >-  \'3in/      -    ^,:cos^  -  I 


(16) 


u 


Substituting;  into   (15) 


L  -  WcoslJ 


(17) 


or 


-/ 


-  Vcos/  (i+  il'Or)  dW 


(18) 


(2)  For  minimizing  the  fuel  used,  1^2,  the  changes  of 
variables  found  useful  are: 
Defining 

Then 


d(   )  S  (   )• 

dx       ■ 


And 


V  -  dV  -  dV  dx  +  V'Vcosl^ 
^      dx  dt 

^  «  d^  -  d^  dx'  -  ^»Vcoslf 

dt       dx  dt 


Substituting  into  (6)  and  (10)  gives: 
L  -  W  v2  ;<•   coslf  -  Wco3j(=  0 

i 

F  -  D.W  V»  Vcos^  -  Wsin}^  »  0 


J 


(19) 


(20) 


(21a) 
(21b) 


The  integral  equation  relating  x  and  W  in  simplest  form  for  minimizing 


•jWg  is; 


As  before 


dV7  -  -  2L£  y 
dx  Vcosft 


(22) 


and 


J       VcosT 


dx 
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Solving  (21b)  for  F  and  substituting  gives 

W  -      I  -  Oi   {Ti  *  WviVcos>(     ♦  Wsin^)dx 

V  cosi$ 


or 


y  '[v^l^     —      —V-J   ^  (23) 

Now,  although  equations  (18)  and  (23)  relate  the  variables 
of  the  orobiem,  the^.  do  not  specify  the  flight  path,  Thej'  are  equations 
for  the  general  motion  of  a  body  in  the  vertical  plane,  expressed  in 
such  a  way  as  to  give  the  relationship  between  range  and  fuel  ex- 
penditure for  any  flight  path.  The  variables  under  the  integrals 
are  dependent  on  the  path,  yet  the  equations  as  such  do  not  specify 
a  path.  To  determine  equations  for  the  particular  path  which  will 
give  the  greatest  range  or  least  fuel  expenditure  \t   is  necessary 
to  look  to  the  mathematics  of  maxima  and  minima*  Presenting  the 
problem  in  more  mathematical  terms,  the  unique  path  desired  is  that 
path  giving  the  extremal  value  of  an  integral.  The  problem  of  finding 
an  estremal  value  belongs  to  that  branch  of  mathematics  known  as  "The 
Calculus  of  Variations",  A  brief  discussion  of  this  subject  follows. 

C)  THE  EULER  EQUATIONS  OF  THE  CAl CUIUS  OF  VARIATIONS 

Let  I_  be  a  quantity  defined  as  an  integral  function  of 
several  variables  as  follows: 

I  =  /a(x,y,2, ,y',z', ,>  dx 

where   (   )•   indicates  the  first  derivative  with' respect  to  x«     The 
Calculus  of  Variations  is  concerned  with  finding  the  extr^nal  values 
of  the  quantity  so  define d^o     The  necessary  and  sufficient  conditions 


4 

FAMILY     OF    PATH  5    FROM 
POINT  "0"  TOT 


1^ 


Tn'W  -  V{x)  +€ni^«(x)  (26) 

With  the  conditions  I 

\(Xo)  -  l(xo)) 

«  f      i.e.    V7(x)=0  at  the  end  points 

Where  (^ia  sl  small  constant  quantitj-  and   T?(x)   is  an  arlitrary 
function  of  x  that  is  continuous  and  has  continuous  first  and  second 
derivativBs  in  the  range  of  interjration. 

^(x)  therefore  defines  all  possible  paths  of  integration 
between  X^jY^  and  X^,I]_  so  that  (2ii)  may  be  written! 

In  -      /  G  [x,Yn(x),VU)]   dx 
or  replacing  Tq     and  Y^'   from  (25)  and  (26) 

^n  "      fGL^>Y(x)  *€n>^(x),Y'(x)  +€j^y^'(x)J  dx  (27) 

X   Y 

It  is^  clearly  obvious  that  I j^ differs  from  Ig  only  through  the  value 
of  £  and  it  may  be  concluded  thatt 

In  -KO  -  f(t)  only 
and  since  €  is  not  a  function  of  x,  the  ordinary  maxima  ahd  minima 
rule  may  be  applied,  ioe<,» 

.S.I   »  0   for  an  extremum 
and,  furthermore,  the  differentiation  may  be  made  under  the  integral,  thus: 


dG  ' 


dx  (28) 


From  (27) 


^        ^G        ♦  d  rY(x)  *    Cn^Cx)] 
de  <)LY(x)  ♦  €^Yy(x)]         d^ 


Since 


dx-dT-dj-dl*   -d7»   -Q 
de       d€       de       d€         d? 


then 


-  •  ^W  * 


60 


or 


d€       d\j{x)  ♦  €n    ij(x)J         '  c)[Y'  (x)   +  £„  ^  (^Oj         * 

<^  "     /    r  c)G  n       ♦  <^G  yp»  T 

^       /    L^(^*^n7^  <^(^  *^n  ^')    J 


.  7?'(x) 


dx  -  0 

(29) 


'Vo'»  'o 


and 


dx  =  0 


From  whence 


^    Y7dX 

Jin' 


^  -O'dx 


Integrating  the  right  side  by  parts  and  remembering -that   >?  (Xq,Iq) 

15 —  ( 


3^,T  "^     '    /      ^^^/    y?  dx 


(30) 


Substituting  (30)  into  (29)   gives 1 

dl  -    I  [dG  ^  «1_.^G  1>?dx  «  0 


rc)G  -  _d ^G  ^y. 

[^       dx    ^j' 


^       /   ["3^       dx    "^' 

\„Y, 

Evaluating  (31)  for  €  «  0  to  find  the  extreraum: 

vX 

.-si"  If  It.]   7*  =  ° 


(31) 
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But  this  is  true  onij  if 


[: 


31   dx 

throughout  the  range  of  integration.  Since  r?  is  conpletelv  arbitrary 
at  every  point  except  the  end  points,  it  necessarily  follows  that  for 
an  extremumt 

jG-Jif^^YO  (32) 

This  equation  is  called  the  Buler  equation.     It  must  be  satisfied  if 
the  integral  is  to  have  an  extreraum  value  and  it  is  through  this 
satisfaction  that  the  path,   Y(X),  is  found. 

Proceeding  to  the  more  general  problem  of  several  variables, 
the  integral  mxist  satisgj   the  definitions  as  before  andi 

I  -     I  Q^,Y,Z, ,   V,Z\ ]  dx 

Proceeding  as  before,  we  arrive  by  an  analagous  deveiopement  at  the 
Euler  equations: 


c^G  -     d  fdO]  -  0 

M     Hx  \M*I  =  0 

3G     -     d    (^G       «  0       etc. 
X7      S"^) 


(33) 
J 


These  equations  (33)  must  be  satisfied  siriultaneously  for  an  extremal 
value  of  the  integral,  such  satisfaction  therefore  defining  the  desired 
path  as  in  the  two  dimensional  case. 

It  should  be  remembered  here  that  an  extremal  value  for  the 
integral  may  be  either  a  maximum  or  a  minimum  value.  Thus  in  deter- 
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n 


mining  a  solution  for  equations  (33) >  we  must  deterirJLne  by  other 
means  whether  a  particular  solution  will  yield  a  maximam  path  or  a 
minimum  path.     For  most  tjpes  of  engineering  work  this  will  be  in- 
tuitivel;s  apparent* 

D)     THE  EULER  EQUATIONS  FOR  FJIXIMUK  RANGE 

a.)  Minimum  Fuel  i-quations   (i.e.,   x  independent  variable) 
Equation  (22)  is  of  the  form 

-  -^2     ■  ^(c)  =  y  Gdx 

•A 
where  1(c)  ■  -W  the  quantity  (of  fuel)  to  be  minimized  and  G  =  CJp  F 

VcosT^ 
the  integrand  function. 


Now: 


G  -  G(Ct,V,i(,F) 
F  »  F(D,V',W,V,^) 
D  =  D(L,V) 


thus 


G  -  G(W,V,7(,0r,^«,V') 
Now,  Cm,  the  thrust  specific  fuel  consumption^  is  a  variable,  and  a  very 

important  variable  where  ran^e  is  concerned,  but  treating  it  as  such 
in  the  derivation  of  i^uier  equations  gives  rise  to  extremely  com- 
plicated relations  and  to  derivatives  on  Cip  that  are  ver}-  difficult 

to  evaluate*  We  have  therefore  elected  to  treat  Gj.  as  a  stepwise 

variable  constant* 

Thus      G  -  G(W,V,7(,V «,(<') 


1? 


And  the  Euier  equations  of  this  function  become: 

1.)     ^G  -  ^d/'c^G  \     '  0 
"3^      .dx  \WJ 

2^)    c^G-     d/^G\    -0 
3.)    dG  -     d  fhG 


iG  -  _d /dG    >|    =  0 


!•)   Derivation  of  the  first  Euler 

c)G     -  _d/c)G  ]  -  0 
"3^       dxV3W 

Now  from  (22) 

G  -  C^  F 


Then 


And 


Vcos^ 

<^G     -       vcos)^  Op  ^    -  Or  F  cos/ 

^  —  —  —       "^  

(Vcos^)2 
-  G  r^F/^V    -  ll 

^G     -_2[ /Af  ^ 

"dV'     vcostj\;^it/ 

d  ^G   =  Vcos^  Or  -f^  ^Xr    -  Oi/^]fVcos?(^  V^?i''sin/) 
^     (Vcos?()2 

Using  an  abbreviated  notation:  sJL  =»  Y.7 

dZ  ^ 

d  fd^Y  C  fFy    -    Fv    (V'cosl^    .  Vi^'sin^l 
dx  V7v7         [jT-        p^^^r^ -j 


or 

ofyi       •      G 


Ffv,'   -  VFy,    +  Fy,  ^'tan^n 

[f —     TT"    r~  J 
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I 


G  Tyt 


rjT"     V" 


Nowt 


GV     -  Gyi  •   "  0  -  Q 


Multiplying  through  by  F  giVes  us  the  first  Euler 
♦  a 

Fy    -  F  -  Fyi*  4  V  Fyi  -  Fyt  ^'tan^  -  0 

2«)  Derivation  of  the  Second  Euler  equation 
(io  -     d    do       ,0 


as  before: 


and 


then 


G  -  Op  F 
Vcosa 


(^G  •     ^^°^^_5(4^ '*'  ^  ^  Vsin/ 
~5o  (VcosV)2 


G 


dG 
17 


F^   -»•  tan/l 
/  J 


d 
Si 


©  -    0,1 


(314) 


Vcos^    ^'    -  F>,/    (V»cosl^    -  V?^«sin^) 


W 


Tvcossys" 


d    <)g 


"T"         "F"\V  / 


G       F,, 


F^.>  -     V«     +  ^»tan/ 
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Now  forming  the  Euler: 


G^  .  G^.'  •  0 


F^  +  tanV 


-  G  F, 


F/  -  V»  +  T(«tan/1 


^i' 


Kultipljing  through  by  F/Q  leaves  the  second  equation* 

F^   ♦  FUnl^  -  F^,»   ^  V»  F^.  -  F^.  V'tanT^  -  0 

3.)  Derivation  of  Third  Euler 

dG  -  _^/^G  '\"  0 
Jw       dxlJwV 


and 


G 


Or  F 
VcosV 


<^G  -     Ct 


Vcos 


'5(\fw/ 


d    f^G\'  0 
dx   UWV 

Qw   -  Gv."'  - 

GfFw 


M-° 


(35) 


(36) 


Therefore  F]^  •  0     and  there  is  no  Third  Euler  in  this  case, 

b)  Maximum  Range  Equations   (i«e,     W  is  independent  variable) 

Notet     d(   )  5  (   )« 
dW 

Frraa  equation  (15) 

G  ■  -  Vcos/ 
Cx  F 

1#)  Derivation  of  the  First  Euler  Equation 


G     -     cos/r»F  »  v('\Z) 
.     G   fl  -  Fv] 
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Vcos^f-    ^F/iV* 


Op  f2 


d    /^G  \  -     Cn,  F^  fv    cos^  Fvi   ♦  VcosT^'Fvi' 

-/vcosJ^Fy,    2F  F«  Ct\ 
V        (^  F2)2  i 


-  VTJ' 


'sin^FviJ 


V»   Fy,    cos/    +     Vcosji'Fy,«  -  Fy,   V  T^'sin/ 


Op  f2 
-  2Vcos/   Fy,   F» 

;3 


Oj,  F2  Orj.  f2 


Qv»' 


Gj  F 

-  G   r  V»  Fy,   ♦     Fyi*     -    -^'tan^    -  2Fvi  F» 
1    V 


F 


J 


G  rFvirr2F»   +  -^'tan^  -  Fyi'     -     V» 


F         F 


'V« 


■] 


Forming  the  Euler 
Gv 


.  Gvi     -    gFJu.  -    Fv  ~  -     G  Fy.      2F«    *    Y»tan^ -JIvlL -J 


And 


F  .  Fy  -  2Fy,   F»    -  Fy,    ll*  Xazi^  l-  Py,  »    +   V»   Fy,    -  0 


(37) 
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2«)     Derivation  of  the  Second  Euler  Squation 
dT  dW   VTF/ 


IT 


V     f-  F  sin/   »  cos  ^  (t)F/^V  ) 


-G 


tan/  +     F/ 


ci  Ci       -     Vcos/    r«    tfF/->/  ' 


TT 


^ 


(F)^ 


-     -VF.i    cos  / 


cr  F^ 

Vcos^F.,*     •  VF.   If'sinji' 


fT 


F2 


V  F^i  cos  V(2F  F') 


d 
dW 


^^  L^    F  F  F  T^2 


Now 


^v  -  °/ 


ftan^rV  FV,  I    +  G  [  V*  F^,    +  jy   -    ^Fy,  tany 
[  F  J  L~  ^^    "^  ^^ 


•  2Fv'FM    -  0 


F2 


From  wtiich: 


F  tan?<  +  Fy  -  V^    F^     -  F^'    +  ^*T  ^  tan/*  2F  ^'F'    -  0 


(38) 


V  F 


Thus  equations  (3li),  (35)^  (36),  (37  )#  and  (38)  represent  extremal 
paths  in  terms  of  flight  path  angle  (5),  velocity  (V),  horizontal 
distance  (x)  and  pounds  of  fuel  (W) •  The  first  three  equations  re- 
present extremal  paths  for  mininum  (or  maximum)  fuel  consumption  for 


*:u 


traveling  a  given  distance  x,  and  the  last  two  equations  represent  an 
Bxtranal  path  for  raaxinium  (or  rainimum)  distance     x    for  a  given  fuel 
expenditure    W» 
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PART    III 
THS  EOUNDSD  CASS 

The  equations  derived  in  section  D  of  Part  II  apply  to  a 
more  or  less  general  propelled  bodj-'  with  unlimited  thrust  and  thus 
unlimited  velocity  and  flight  path  angle.     They  do  not  consider  any 
boundaries  on  the  tiirust  available  to  follow  the  extremal  path  that 
they  designate.     This  ext-remal  path  maj  very  i^erll  call  for  altitudes, 
velocities,  and  accelerations  above  the  physical  ability  of  a  parti- 
cular airframe  -  engine  combination  for  ail  or  part  of  the  path.       ** 
Tho&e     equations  are,  therefore,  the  absolute  extremums  for  an  ideal 
body  without  limitations  imposed  on  any  of  its  parameters.     They 
define  a  path  that  is  the  best  possible  and  should  be  followed  when- 
ever that  path  is  within  the  operating  capabilites  of  the  particular 
propelled  lody. 

The  primarj  limitation  imposed  by  physical  factors  in  any 
particular  case  is  that  of  thrust.     Thmist  varies  from  zero  to  a 
fin-ite  maximum  value  that  is  a  function  of  the  individual  thrust 
producer  and  other  variables.     For  a  given  ".type  power  plant  the 
maximum  effective  thrust  is  usually  a  function  of  the  variables  of 
the  problem©     For  example:  The  thrust  of  a  rocket  is  a  constant^  the 
thrust  of  a  ram  jet  varies  with  velocity  and  altitude;  the  thrust  of  a 
turbojet  varies  with  velocity  and  altitude.     Therefore,  for  any  ^ven 
airframe  -  engine  comlination  there  is  a  limiting  envelope  or  boundary 
of  thrust  which  limits  the  path  variables  such  as  velocity,  flif^ht 
path  angle,  altitude,  etc.     V/henever  the  general  path  equations  of 


I 


I 

II 
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Part  II  demand  a  thrust  that  exceeds  this  boundary,  the  path  becomes 
unrealistic  and  of  academic  interest  only.     It  then  becomes  necessary 
to  compute  a  new  path  from  this  point  that  will  be  an  extremal  based 
upon  a  maximum  continuous  value  of  thrust  available.     Euler  eqiations 
will  now  be  derived  to  define  an  extremal  path  for  the  case  of  b9unded 
tlirust. 

A.)     The  Bounded  Case  for  Kinimura  Fuel  for  Turbo  jets 

As  was  said  before,   the  thrust  of  a  tiirbojet  varies  with 

altitude  and  velocity,  i.e.: 

F     »    f (V,h)     where  f  pertains  to  bounded  thrust 

But  since  altitude  is  a  function  of  path  angle,  we  may  also 

say: 

F    -   fiv,i) 

From  before: 

G  -  2Llj 

Vcoso 
But  now 

G     -     Or  f 

Vcos^  (39) 

This  function  is  ind^endent  of  the  characteristics  of  the  airframe  - 

engine  combination  and  must  be  modified  to  include  them©     From  (21b): 

V  cos^    »     g(f  -  D  -  Wsin^) 

Substituting  this  into   (39)  above 

G    -     Or  f  W  V  (1;0) 

g(f.D-Wsin^} 
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!•)  Derivation  of  the  First  Buler  Equation 

<^Q  -d.  /<^a^  -  0 


G_  -     Pr  f  w -  _G_ 

.V»     g(f  -  D  -  Wsin^  )       VI 


But 


Then 


Now 


W 


^G   -   W 


d  /  jiG  \  -  yiw^  ,  W»  V",  (I42) 

dx  l"SF;       (v»)2 

Forming  the  Euler  (ia)  -  ih2)   -  0 


/_iG  '\^  ■  W  ffv  -  fV  -  ^V   1- 


aG  -  d  /  ^Q]     -  W  I  f V  -  f V  »  Dy    |,  ytw^'  ,  w'V"  -  0 
"W       dx  l"3Vr;       \^       (fJLWsinV)!      (v«)2 


From  which 

-0  (^) 

J 


_W^  -  V  ("fv  -  f Y  »  Dy   -1,  yw 
W      [f     (f-D-Wsin^)J   "f 


2.)  Derivation  of  Second  Kuler  Bquation 
^)^       1.  f       (f  -  D  -  WsinV  ; 


ihh) 


dV  (f-D-Wsini) 


J    (^h§\  -    W»  r(f-P-V-^sin^)D^'     -    D^   (r-D'-W'sin^  »  Wy'cos^)"] 
dK\^i"j  [ifliSIwiiir^  (f-D-WsinT()2  J 

*  ^'r_£il_l  (1*5) 


11 
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Fomin-  the  Euler   (kh)  -  (16)  -  0 


■■[? 


*   w 


=  0 


From  v;hi  ch 


TT  f 


Partial  derivatives  of  drag  with  respect  to  V   and  fnrust  with 
respect  to  li'   are  zero  since  drat^  and  thrust  are  not  functions  of 
V*   an-^  ii*   resp active! J-. 

E)     The  Bounded  Case  for  Kaxinum  Range     (i.e.  with  W  as  independent 
variable ) 

As  before 

G     =»     Vcost( 


Gp  F 


and 


so 


F     -     f(V,h)     -     f(V,^)     only 
G     ■     Vcosft 


(1:7) 


Introducin'^  the  airframe  characteristeics: 
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From  (13a) 


L*WVi^'CtF     -  Vcosi     -  0 

g 

cosl^    -     L     +     V  1^»   Or  f 


g 

Substituting  this  into  (li?)  leaves i 

OrWf         g 


.-It 


1«)   Derivation  of  First  Euler  Equation 


^G-         L         ♦VLv.VLyf*       2VY'  {U?) 

c>G       -    _Xil'      (thrust  is  not  a  function  of  V) 


L 

-     V  Lyf 
GpWf2 

"♦• 

2  vy 

C^Wf 

g 

-1-  /"  ^G  ^     ■     ^^''  "  ^'^^'  -     V  Lvi(Wf'    ^  f) 

dv;  V"^/  Or  w  f  Or(wf)i^ 


Forming  the  Euler  (1^9)  -   (50)  »  0 


(50) 


L  +  VLy  -  VL  f V  +  2CjW7j»fV  -  VL^i »  -  V'Ly,    +  VLvtJwf'    +  f  ]  °  0 


2»)  Second  Euler  liquation 

£]  ^*  '     (52) 


ha     - 

OrWf 

J2-  - 

Gj.Wf            g 

II 
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(53) 
Fomdng  the  Euler  (^2)   -  {S3)  "0 


L|  -  L  f^   -  L^»   -  V   L,'    +  Ly'[wf«    +  f ]  -  20rW  f  V«      -     0 
~f"  ?"  Wf  "g" 


ish) 
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PART  IV:  EVALUATION  OF  THE  DERIVATIVES 
AND  FORMULATION  OF  THE  FINAL  EQUATIONS 

The  Euler  equations  as  developed  in  Parts  II  and  III  contain 
derivatives  on  several  of  the  variables.  To  determine  the  desired 
path,  these  derivatives  must  be  evaluated  and  included  or  discarded 
depending  upon  their  effect.  In  this  thesis  only  the  derivatives 
for  the  minimxDn  fuel  path,  both  bounded  and  unbounded,  will  be  deter» 
ndned*  The  evaluation  of  the  derivatives  for  the  maxiinuin  range  path 
(W  independant  variable)  is  analogous  but  the  expressions  are  much  more 
unwieldy. 

Cwisidering  first  the  derivatives  to  be  evaixiated  for  the  first 
Euler  in  the  unbounded  case,  we  havej  (l)Fy,  (2)Fyi,  {3)Tjjt**     From 
equation  (21b) 

F  -  D  +  W  V«  Vcosl^  -^  Wsin^  {$$) 

i 

Taking  the  desired  derivativeaj 

Fy  -  Dy  ♦     W  V»    cos^  (56) 

Fyi-  Dy,  *    W  Vcos^  (57) 

Fyi »   »  Dy, »   "»•  _WM[cos^     +  J^V»cos^  -    i*  W  VsinO 


g 


•  D^,»   +    W  Vcos/Tw   +     V»  -    •^•tan^l 

T       Lw       V  J 


(58) 


Consider  now  the  derivatives  to  be  evaluated  for  the  Euler  in  jr 
for  the  unbounded  casej  the  following  appear:  (l)F^,  (2)F^< , 

(3)F^i\  From  (55) 

F  -  D  •»•  W  V  V'cos^  *  Wsin^ 
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Whence : 

^i  '  ^i  "^  JL  ^  V'sin^  +  Wcos -i  (59) 

g 

F^.  -  D^i  (60) 

F^.  •  -  D^.«  '  (61) 

These,  then, are  the  desired  derivatives  as  functions  of  drag  or  of 
drag  derivatives.  The  drag  derivatives  will  be  determined  after  a 
look  at  the  derivatives  appearing  in  the  bounded  equationg.  From 
equations  {h3)   and  (U6)  for  the  bounded  case  we  have  the  following 
derivatives: 

(1)  W* 

(2)  W" 

(3)  fv 

(5)  f  • 

Now  W  -  Or  f  W  V* 

g^f-D-WsinlO  '   (^2) 


From  which 
W"  » 


Or  f  W  V»  r  f  I  +  w'  +  V"  -  f^'  -  D«  -  W'sin'J^  -  W)i'«cosY)] 
jCf-DLWsin^)[T"  "TT  "TT       (f-D-WsinS;        J 


(63) 

Since  the  thrust,  f,  used  in  the  bounded  equations  is  an  actual 
thrust,  (i.e.  thrust  available),  engine  data  for  the  specific 
engine  -  airframe  combination  will,  supply  the  information  necessary 
for  evaluating  the  derivatives  on  thrust. 

f  ■  f(V,h)   A  known  relationship  for  a  given  aircraft 
Then 

fy  »  ^f/cIV  This  can  be  taken  directly  from  engine  curves  (6l4a) 
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.^i 


Now 


f,  -  ^f   ^h 


(6l4b) 


and  ^f/^h  is  available  from  engine  curves  while  hh^)i   is  a  simple 
geometric  relationship© 

f »  -  ^  V»  ♦  _^  V  (6Uc) 

To  evaluate  the  derivatives  on  drag  it  is  most  convenient. to 
express  drag  in  the  following  forraj 
D  -  D^  +  D^ 

•  _i-/^S  v2Cd^  +  1  /?S'#  q^ 
2  /  2  /     "^TeJir 

Or,  substituting  for  C^  its  equivalent  in  terms  of  lift  from  equation 
(21a): 


D-    1    /)S  V^Cd    »  2   Odcosi)^    /l-*-    V^Y   y       ^ 


"T" 


/^ 


ireAR/)sl^ 


(65) 


:/>; 


In  taking  derivatives  on  the  above  equation  for  drag,  the  subsonic 

region  only  will  be  considered.     No  Mach  effects  will  be  present 

then  and  both     Cn      and    e    will  be  treated  as.  constants.     An  extension 
^o 

of  this  analj'sis  to  include  Mach  effects  will  be  found  in  Appendix  B, 
Thus 


'  ireAR  z>S  V3 


r- 


=     2  Do  -     2  Di 


1  - 


(66) 


1  + 


8 


V    -0 

Dv»*   ■  0 


(67) 
(68) 


I 


31^ 


^i      ■    -h  W  cos)(  sin//l  ♦     v2  5'V 


^5f' 


-2  D^  tan^ 
U  (Wcost)2  f 


g 


(69) 


Dy,. 


-      2  Di 

W2 


,       V2^ 


g 


(70) 


cos^^    (l  >  1!^J  2w,   .  2^'tan;(  *  {7^"  *  2P 


ireiR  /)S  g 


:V^K"  •«'  2/^'VV*)--    /)t 


2Di 


2W«   -  2  j'tan?-  -  £J^  +  V  (V^   *  2V'jr') 


"VT 


1  *  Pi^Sv 


(71) 


2DqJ"  y   ♦    /?«]    -  2Dj^ 


V         ^J 


g 


XI 


-      W    -^      iD« 


S 


,  /O  ,  and  its  derivative,  /O  »,  may 


(72) 


The  atmosphere  parameter 
be  evaluated  as  follows} 

For  altitudes  below  36,089  feet: 
pAl  -  6.89  X  10-^  h  )^^^^ 


r-i 


t 


'   -    -I1.256  /)o(l  -  6.89  X  10-^  h  )3*2^6  (5^89  X  10-6     ^ 


r 


"ar, 


But      dh     -  tanV 


Then: 


/>»     -    -  29.38  tan^ 
Y>  10^  -  6.89h 


(73*) 


And 


For  altitudes  above  36,  039  feat: 


r-  r 


X 0.2971)  e^^^^^ 


rbc 

^oTHT 


7^ 


-  tanV 


(7h) 


Tbf!  formulation  of  the  final  equations  for  thf»  unbounded 
caso  nay  no-'f  be  made  as  follows: 

'   -  0 


Fv  -   F 


T" 


V 


Substutin-  for  ?y,   F/^/,  F^i',  and  F^i   fron  (55),    (56),    (57), 
and  (53);  and  for  Dy  which  occurs  in  Fy^  from  (66)  gives  after  rearranf^ing: 
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^o  -    Pii   (^  "••^^i   U    ^'    Vcog)J' 


ill  *  ( _!.  ■*■  2^*1    *-"ii  i 


-  0 


'i^ 


(75) 


rron  (35),   tlie  second  Eulsr  equations: 

Fy  +  FtanV  -  F^\»    -  P^l(»tani  -    V*   "j     «  0 

^ubstitutinfT  for  F^,  Ftant,  F^(»   and  F^«  from  (59),    (60)  and  (6l) 
and  for  the  dra?^  derivatives  which  aro  introduced  from  (69),    (70) 
and  (71)   fjives  after  rearranf'in^: 


(Do  -  Di)   tan^  +     \l 
cosjf 


-     7»    +     V     ^VV"   +  27«r'  ) 


2  TTj  j  2  y 

\   n 


-.^»tan);.£» 


(76) 


Solving  (75)  for  V   and  (76)  forV'  yields  the  two  equati 
whose  solution  will  define  the  optimum  path. 


ona 


WOrv2co8^    )     \        ~ 


t    5^ 


2Di 


♦  gOp  •».  2gl(» 


tano 


J   ^ 


(Do  -  Di)  tani  ♦    W 

cos? 


I    V2  i 


(77) 

r 

2W« 

"V- 


-     V«  -  ^'tanl^.  r« 


^. 


(78) 


The  final  equations  for  the  bounded  case  are  formulated  in 
like  manner: 

From  (li3)»  the  first  JEuler  equation 


W"     •  V»  i     ^V    •        f  V  -  Dy 


i     -V     •         *v V 


-     V»'     -     0 

"IF 


Substituting  for  W"A'*  from  (62)  and  (63)  and  for  Dy  and  D«  which 
then  remain  from  (66)  and  (72)  yields  on  rearranging: 

+    W»   -  _^  I    -  f  •  ♦  W'sinft 


2  D, 


V 


1     *    />* 
V        1^ 


'-^'•tan^  ♦ ± 


W^i 


T    KTj 


+  W'Vcos)i    +    /^  -  ZI£v  "^  Jil^     (f  -  D  ^Wsin)^)  -  0       (79) 
(^   f         f  w    y 

From  (U6),   the  second  Suler  equation^ 


-b[l  D^*  *  -£l  (f  -  D  -  Wsin^O  •  f ^  -     D^  -  \lcosi  *  D^i» 
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T^  .    (ft  .  pt  »\i*Bi.Tii  ^WI^'cobV)     -  0 
^  (f  -  D  -  WsinTj;  ^ 


Substituting  for  W'/W*   from  (62)  and  (63)f  for  ly  ,  Djji%  and  D» 
from  (70),    (71) >  and  (72)   jlelds  on  rearranging j 

f»   +  3W«   ♦     V"  -     2(f»  -  W»sin^  -  W  ^li'cosK)  -  2D» 


(f  -  D  -  Wsin^) 


— ? 

♦     V(V^"  -»•  2V'  i')     -    Pj^  -  2  T(»tan^     -  0 


M    (f  -  D.Wsinli)  -  f  )^  ♦  2Di  tan^  -  WcosX 


77  ? 


(80) 


Solving  (79)  for    ^"  and  (80)  for  V"  yields  two  equations 
whose  solution  will  define  the  desired  path J 

^'"    ■  i\     1    1  ff'  -  W«sin^- W  i^cosi  -    /  f»  .     V'fv    +    w»  )  (f-D-WsinT() 

[T-    -r    ^yj      v-ir    -r/   /i^  6'vg\ 


r 


♦    i'^tani^J^*Jl){  ^*V) 


(81) 


V« 


'-  V«f»  -  3W'V'   +     2V»(f»   -W'sin^.  -  W  ^»cos^)  -  2V'D» 
"r~      "^ (f  -  D  -  Wsin^) 

-        ("^  )       \   ^i   (f  -  D  -  Wsiid)  -  f^  +  2Ditan^  -  Wcos/j 


Wi    {Vi"*  2V«  VO     +     VV'»     +  2V'  )l'tan]( 


(82) 
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PART  Vi     METHOD  OF  SOLUTION  OF  TH2  SQUATIQNS  FOR 
AN  EXTRH^IAL  PATH 
Considering  first  the  unbounded  or  general  equations  for  an 
aircraft  fijing  a  minimum  fuel  path,  the  following  equations  control 
the  aircraft's  motioni    (21b),    (65),    (77)  and  (78).     These  equations 
will  be  repeated  here  for  ready  reference: 

A,  The  general  equation  of  motion  along  the  flight  pathx 

F  -  D  •••    W     V«   Vcosi  ♦  Wsin^ 

T 

B,  The  drag  equation: 

D  -  1  /7s  V^Cjj^   ■».  2(Wcos^)^  /  1  +  ^i*Y 
"^/      *    1fe«l/»S  i^        \  g  / 

C,  The  first  Euler  equation: 


(A) 


(E) 


A 


v 


A 


Dofl-  VCt\  .  1^ 


WOpV^cos^    /        \  g 

g  Y  *  g^  +  2g^^» 
D,  The  second  Euler  equation: 


tan)i 

J  -w 


j 


(c) 


-^" 


("o  -  %)  ^"^  *  _S 


-  /_i_  *  ^'1  ^ 


where 


-    V«  -     i'taii^  -    PI 


-     2V»'i(' 


v2     yL 

(D) 


W 


^0-4-/ 


1    pS  V'^  CD 


*Dj^  -     2(Wcos'?i')2      /l  +  jW2 
remPs  l'^     \  g     ^ 


RP 


This  system  is  composed  of  four  nonlinear  differential  equations 


liO 


in  5  and  V«  Three  of  the  equations  are  first  order  and  one  is  second 
order.  This  is  equivalent  to  a  third  order  equation  in  one  unknown 
and  thus  three  initial  or  boundary  conditions  may  be  specified. 

Since  no  general  closed  form  solution  exists  to  this  system,  we 
are  forced  to  integrate  by  means  of  a  step-by-step  process.  This 
requires  the  use  of  initial  conditions  rather  than  boundary  conditions. 
From  equations  (A)  and  (E)  it  can  be  seen  that  the  forces  on  the  air- 
plane depend  only  on  V,  V* ,  ^  ,   and  t '  and  therefore  our  initial 
values  should  be  given  in  terms  of  three  of  these  at  a  given  altitude,  h« 
Since  V  can  also  be  determined  directly  in  terms  of  the  other  three 
from  equation  (C),  it  will  be  most  convenient  to  specify  V,  o  ,  and 
i^   initially. 

For  a  given  aircraft  we  know  Wq  and  Cp  ,  and  for  the  unbounded 

case  we  assume  unlimited  thrust  and  knowledge  of  Cj  at  any  thrust. 
Now,  if  we  specify  V,  "jf  ,  and  3»  initially,  V*  is  computed  from  (C) 
and  D,  D©  and  D^  from  (B).  With  V*  and  D  known,  F  is  determined 
from  (A).  Since  Cj  is  a  function  of  F^  it  was  necessary  to  estimate 

initially  to  find  V,  Now  that  F  has  been  determined,  an  iterative 
procedure  should  be  used  to  correct  V*,  Gj.  and  F.  W  is  found  easily 
from:  , 

W»  -  »  Cy  g  (E) 

Vcos'lf 

Now  1^"  can  be  ccwputed  from  equation  (D),  Thus  we  have  defined  all 
quantities  for  our  starting  point  by  initially  specifying  three 

variables  V,  a  ,  and  o  *o^ 


hi 


By  using  the  first  two  terms  of  a  simple  Taylor  series  ex- 
pansion we  may  now  integrate  step-wise  to  point  number  1  and  so  on  in 
the  following  manner x 

Knowing  Vq',    %  q* ,  Wq*   and  i  q"  we  compute  by  Taylor  series 

^1     -     Vo  >  Vo'Ax 

^1     '    io  *  ^o'Ax  ^   ^ 

.  ^  (F) 

Wi     -    Wo  +  Wo'  Ax 

Also  since 

h»     -     tan^ 

h;L    -    hg  ♦  h'  Ax 


(F) 


Now  with  V;l,    iif    ^i,*  f  W^,  h]^  and  equations   (A)  through  (E)  we  can 
determine  ail  quantities  at  point   (l)«     Using  Taylor  series  again  we 
proceed  to  point  (2)  and  so  nno 

If  greater  accuracy  is  required  than  obtained  by  two  term 
Taylor  series,  additional  terras  of  the  series  may  be  used,  the  Milne  method 
may  be  used,  or  the  increment  in    x    may  be  decreased.     Decreasing 
the  increment  in    x    reduces  the  errors  due  to  non-linearities  of 
the  function  involved.     The  Milne  method  is  described  in  reference  2 

In  performing  a  step-bj-step  solution  to  the  unbounded  case, 
several  considerations  should  be  kept  in  mind  to  insure  that  the  path 
so  determined  will  be  a  "flyable"  onet 

1,)  The  equations  can  possibly  give  a  maximum  fuel  path  instead 
of  a  minimum© 

2,)  Flight  path  angle,  a, cannot  be  allowed  to  drive  altitude  to 
negative  values  or  the  a/c  into  stall. 
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3»)  Accelerations  and  decelerations  must  be  within  the  physical 
capaLilities  of  the  aircraft.. 

h»)   Minor  oscillations  of  the  computed  path  in  its  initial 
stages  are  probably  due  to  computational  variables  such  as  the 
size  of  X  and  to  the  initial  conditions,  and  can  be  faired  out  for 
an  actual  flight  path©. 

Juggling  of  initial  conditions  and  increment  in  x  will  help  avoid 
these  difficulties  at  the  beginning  of  the  path* 

To  compute  a  minimum  fuel  path  for  the  bounded  case  equations 
(C)  and  (D)  must  be  replaced  by  (6l)  and  (82).  Equations  (8l)  and 
(82)  will  be  restated  here  only  in  functional  formt 

First  Suler  Equation 

f(V,  i  ,   V,  ^S  i")   -  0  (G) 

Second  Euler  Equation 

f(V,  i   ,  V«,  )(',  V",  ^«)  '  0  (H) 

It  must  be  remembered  in  the  bounded  case  that  thrust  is 
fixed  as  a  known  function  of  V  and  "i  for  a  particular  aircraft. 
"With  this  in  mind,  a  close  imspection  of  all  of  the  equations  for  the 
bounded  case  (i.e.  (A),  (B),  (G),  (H))  tells  us  that  the  second  Suler 
equation  (H)  is  not  necessary  to  the  solution  of  the  system.  This 
is  as  we  would  suspect  from  the  physical  conditions  of  the  bounded 
problem,  since,  by  setting  thrust  at  a  fixed  rating,  say  maximum 
continuous,  it  becomes  a  known  function  of  V  and  a.  Thus  one  less 
equation  is  required  for  solving  the  sustera« 
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This  can  be  seen  more  clearly  from  a  look  at  (G)  and  (H), 
with  the  balance  equations.  We  have  two  second  order  differential 
equations  with  second  order  terms  on  JT  in  both  (G)  and  (H) 
equations  and  second  order  terms  on  V  in  (G)  only,   jr"  could  be 
eliminated  between  (G)  and  (H)  and  our  system  would  be  equivalent  to 
a  third  order  equation  in  one  unknown.  Thus  three  initial  conditions 
nav  be  specified.  Looking  at  (G)  alone  with  the  balance  equations 
we  also  have  an  equivalent  third  order  equation  in  one  vaMable» 
Thus  three  initial  conditions  may  be  specified  in  either  case. 

If  we  initially  ^ecify  V,  y,  and  ^r*  as  before,  for  a  particular 
aircraft,  F,  Ct*  and  D  are  determined.  Then  since  F  and  D  are  known, 
V'  is  determined  from  equation  (A)o  If  V*  is  determined  from 
equation  (A),  the  V"  term  in  equation  (H)  is  redundant.  The  other 
variables  in  equation  (G)  are  functions  of  known  or  specified 
variables  so  V,  JT,  V*,  5"'  and  jf"  are  related  uniquely  in  (G)  and 
therefore  (H)  can  be  disregarded,  Ey  using  only  the  first  Buler 
equation  all  quantities  are  specified  or  determined  for  the  initial 
point  of  the  boxinded  problem,  Py  equations  (?*)  representing  the 
first  two  terras  of  the  Taylor  series,  we  can  again  integrate  stepwise 
to  point  1  and  so  on  as  before, 

,  The  bounded  problem  and  also  the  more  general  unbounded  problon 
each  have  specific  applications  which  fit  ccxnpletely  within  the  scope 
of  one  or  the  other  of  the  two  situations.  It  is  felt,  however,  in 
the  absense  of  any  computed  results  that  the  general  solution  will 

Ok 

be  composite  of  the  two.  Due  to  the  similarity  of  the  systems  of 
equations  for  the  two  cases,  transition  from  one  case  to  the  other  is 
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an  easy  matter  at  any  point  oh  the  path.     When  thrust  required  ex^ 
ceeds  that  available  while  using  the  uniounded  equations,  transition 
to  the  bounded  system  of  equations  should  be  made.     This  transition 
should  be  accomplished  at  the  last  point  on  the  path  at  which  the 
thrust  requirement  is  equal  to  the  thrust  available* 

Now,  while  proceeding  step-by-step  along  the  path  computed  by 
the  bounded  system  of  equations,  it  is  necessary  that  each  path  point 

/ 

be  checked  in  the  unbounded  equations.  At  that  point  on  the  path 
where  the  unbounded  thrust  requirement  comes  within  the  envelope  of 
thrust,  available  from  the  power  plant,  transfer  back  to  the  unbounded 
system  of  equations  should  be  made  and  coinputations  continued  until 
such  time  as  the  thrust  requirement  again  forces  transition  to  the 
bounded  case.  This  procedure  is  continued  until  the  required  distance, 
X  is  coveredo 


IS 


COTiaUSIONS 

It  is  concluded  that  a  system  of  equations  for  Tnaximum  range 
or  minimum  fuel  flight  paths  ma^  le  written  which  are  not  restricted  to 
level  flight  or  to  negli^ille  accelerations.     The  equations  can  be 
;-^rtnlied  to  flight  paths  which  coirjr-ionce  at  take-off  or  at  anj--  other  set 
of  initial  conditions  from  which  an  extremal  path  is  desiredo 

It  is  further  concludefi  that  bj   specif jdn^  three  initial 
conditions,  nethods  of  solution  bi   stepwise  integration  exist  for 
the  system  of  equations. 

It  is  also  concluded  that,  should  the  extremal  path  call  for 
greater  thrust  than  is  available  from  tlie  aircraft,  a  maximum  range  or 
mininum  fuel  path  can  be  computed  based  on  maximum  continuous  thrust 
available.     There  are  specific  applications  for  both  types  of  paths 
individually  as  well  as  for  a  corposite  patn. 

From  Appendix  A  it  is  concluded  that  classical  parameters  for  > 
maximiOTi  range  are  restrictive  applications  of  the  more  general  equations 
derived  in  this  thesis© 


APPB-imX     A 


Reduction  of  Euler  Equations  to  Classical  Range  Parameters 


1|6 


Frcsa  the  Euler  equations  representing  a  inaxiraum  range  path  at 
variable  altitude  it  is  a  simple  matter  to  derive  equations  or 
parameters  for  maximum  range  in  level  flight.     Starting  with  equation  (77)* 

.     Gj  A   ^  V^iL 'SI  -    Wtan2(/i  >  Ct   V> 


g    \  g 


V2 


Imposing  the  conditions  for  level  flight  (  "/  ■  V' '  ■  O)  we  have 


remaining: 
V» 


2   V     g  /   Ap\y  (    V  .    g 

Substituting  f or  W^- its  level  flight  equivalent  (W  - -^/^sV^Cl) 

2 


_  91^12  -     losV    Q^       -i/^sVCo^ 

TTefiRs/aV^  2'         we^E  2' 


r 


Then 


And 


V'WDr 


osCDo  a  -  Pr  v^  +  ipsvcD.  / -  3  -  Pr 

i^i    ~g-;  2/      [  1  — 


v'Cl  Or  v2    -    odq/i  -  Ct  vV  Cd^  a  *  Ct  V 


g'- 

Assuming  nominal  values  of  Oj,  V,  and  C^,  sa^: 
Op  »  1  lb/lb  hr  -  1/3600  sec 
V   -  liOO  ft/sec 

Cl  -  .2 
»  in  slightly  different  form 


g 


(A.l) 


(A.2) 


m 


Then  Oj  V  •  •0035  and  may  be  considered  negligible  with  respect  to  1» 

e 

And  V^Cj^CxV^  »  »009  V  where  V»  is  in  ft/sec  and  is  an  extremely 

— F~  ~^^ 

small  quantity*     This  tcnn  may  be  considered  to  be  zero*     V/ith  these 
assumptions  the  equation  reducec  to: 

CDq     -     3  CD^     -     3  Cd     -     3  Cr,Q 

Or: 

%     '     k%o  (A-3) 

Equation  (A#i)  represents  a  parameter  defining  the  speed  schedule  for 
maximum  range  at  any  one  given  altitude*  This  parameter  is  perhaps  not 
as  familiar  as  the  classical 

parameter  but  it  can  be  shown  that  the  two  are  equivalent* 
Differentiating  ^^/q^  "   K 

_L-  c£""  -  K  dCD  -  K  ^PL 
2  dpL      "^effl 


Then 

K.  - 

teff? 

and 

CD 

1 

-  .'item 

from  which 

Cd 

'     h  Gl2 

teAR 

and 

-    ^  % 

U  Cd^  -  li  Cd  -  hCDo 


(A.$) 


APPSNDIX     B 
Consideration  of  Itech  Number  Effects 

In  the  foregoing  developement  of  maxinium  range  and  minimum 
fuel  equations,   onlj^^  the  sub-sonic  speed  region  was  considered.     In 
this  section,  Kach  effects  on  the  equations  will  be  determined  and 
applied,  thereby  removing  any  speed  limitations* 

Mach  number  appears  explicitly  in  the  drag  terras  of  both  the 
unbounded  and  bounded  cases,  and  may  possibly  appear  in  the  thrust,  f, 
of  the  bounded  case.     Therefore,  by  determining  the  Mach  corrections 
to  the  D  and  f  terras  and  their  derivatives,  the  canplete  system  of 
equations  is  corrected  and  is  applicable  to  the  entire  speed  range. 

Considering  drag  first,  the  ^5ach  effects  thereon  may  be  accounted 
for  as  follows: 
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Now  Cj)     and  e  are  functions   of  Mach  number.  Defining 
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Taking  derivatives  on  dragt 
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Since  neither  drag  nor  Mach  number  are  functions  of  V,  as  before: 
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0  ,  also  since   ^M   ■   c)^' 
the  derivatives  on  drag  with  respect  to  o    and  a'  are  unchanged 


from  their  subsonic  form  provided  Cj)  and  e  are  inserted  as  proper 

functions  of  Mach  number* 

To  evaluate  Dy',  from  (70)  and  the  above  discussion: 
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To  evaluate  D» 
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As  was  previously  noted,  the  bounded  thrust,  f ,  may  possibly 
be  a  function  of  Kach  number.  If  this  is  the  case,  derivatives 
may  be  taken  in  the  normal  manner  as  follows t 

f  -  f  (V,  h,  K) 
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And 
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The  partial  derivatives  with  respect  to  the  various  variables  being 
obtained  from  known  power  plant  characteristicso 
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